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Abstract 



For a long time, the well-known supercritically accreting binary SS433 is 
being proposed as a prototype for a class of hypothetical bright X-ray sources 
that may be identified with the so-called Ultraluminous X-ray sources (ULXs) 
in nearby galaxies or at least with part of them. Like SS433, these objects 
should be associated with optical nebulae, powered by both radiation of the 
central source and its wind or jet activity. Indeed, around many ULXs, 
bright optical nebulae (ULX Nebulae, ULXNe) are found. Here, we use 
SS433 as a prototype for the power source creating the nebulae around ULXs. 
Though many factors are important such as the structure of the host star- 
forming region and the possible supernova remnant formed together with the 
accreting compact object, we show that most of the properties of ULXNe may 



be explained by an SS433-like system evolving for up to about one million 
years in a constant density environment. The basic stages of evolution of a 
ULX Nebula include a non-spherical HH-region with a central cavity created 
by non-radiative shock waves, an elongated or bipolar shock-powered nebula 
created by jet activity and a large-scale quasi-spherical bubble. 

Keywords: ISM: jets and outflows, ISM: bubbles, stars: individual 
(SS433), X-rays: individual (ULXs) 



1. Introduction 

The issue of supercritical accretion is mor e than three decades o ld. It 



was first considered in the seminal paper by IShakura and Sunyaevl (119731 ) 
who pointed out the significance of Eddington limit for accretion processes. 
A mass accretion rate exceeding the critical, or Eddington, value leads to 
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Eddington luminosity limit violation in the fra mework of standard disc ac 



cretion. The disc may become either advective ( lAbramowicz et all Il980f ) 



or 



outflow-dominated (see iPoutanen et al.l (120071 ) and references therein) if the 



critical accretion rate is exceeded. 

Supercritical accretors are expected to appear among massive black-hole 
X-ray binaries dur i ng th ermal and sometimes nuclear-timescale mass transfer 



( IRappaport et all 120051 ) . The released power is in excess of ~ 10 ergs 



(Eddington luminosity for a conventional 1OM black hole), that makes su- 
percritical accretors in nearby galaxies prospective targets for X-ray, UV and 
optical observations. In different models, however, the exact observational 
properties differ. Both outflows and advection (advective discs are geomet- 
rically thick) taken apart make sources anisotropic and soften their spectral 
energy distributions (SED), thereby increasing the contribution of extreme 
ultraviolet radiation in the overall energy output. That make s supercriti- 
cal accretors possible en ergy sources for bright optical nebulae (JKatzl . Il986 : 



Abolmasov et all 120091 ) 



SS4 33 is known as a unique Galactic object in the supercritical accretion 
regime (IFabrikaJ . J2004J ) . Although its X-ray luminosity is much lower than the 
expected Eddington luminosity (~ 10 36 -rl0 37 ergs _1 , depending on the inter- 
stellar absorption value), similar objects observed at low inclinations < 20° 
are predicted to be objects of high (> 10 39 ergs _1 ) X-ray luminosity (IKatz . 



19861 : IPoutanen et all 120071 ) . It allows to link SS433 and the still rather dim 



concept of supercritical accretion with extragalactic Ultraluminous X-ray 
sources (ULXs) known to have X-ray luminosities 10 39 -i-10 40 ergs^ 1 . Though 
not unique in this sense, the supercritical accretion model qualitatively ex- 
plains most of the observational pro perties of ULXs, in c luding their envi 



ronmental and statistical properties (ISwartz et al 
2007bf ). 
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This interpretation is indirectl y supported by the fact tha t some ULXs are 
identified with powerful nebulae (jPakull and Mirionil . 120031 ) with sizes from 
about 10 to hund reds of parsecs. The optic al properties of these nebulae are 



very diverse (see I Abolmasov et al.l (J2007bl ) for review). Some appear large- 



scale shock-powered bubbles. However, in several cases an Extreme Ultravio- 
let (EUV) ionisation source is needed to explain the properties of the nebula 
or some part of it or some particular emission line s. In the case of the large 
nebu lar complex MH9/10/11 around HoIX X-l (I Abolmasov and Moiseevl . 
20081 ). two regions powered by shocks and EUV radiation are clearly distin- 
guished by kinematics and emission-line spectrum. 



SS433 is surrounded by th e radionebula W50 (IDubner et allll998l) k nown 
to harbour optical filame nts (Zealev et all Il980t Ivan den BerghJ . ll980T) . Re- 
cent spectroscopic results (jBoumis et al.l . l2007t lAbolmasov et all |2010| ) show 
that bright [N II] and [O III] emissions are broadened by ~ 50kms _1 . Their 
relatively high intensities (in comparison with Balmer emissions) are easy to 
explain if Galactic abundance gradients are taken into account. The struc- 
ture of the nebula is seen in details, but it is difficult to infer its integral 
properties such as total line luminosities and line ratios because of heavy 
and patchy absorption toward the nebula and its large angular size (about 
1 deg). The total emission line spectrum of W50 may differ significantly from 
the spectra of the bright filaments. 



T here were several works con cerning the evolution of W50 (IZavala et al. 



20081 ; IVelazquez and Ragal . |2000|) . They agree about the leading role of the 



relativistic jets of SS433 in forming the peculiar "seashell" shape of the neb- 
ula, but its central round core is usually attributed to the pre-existing su- 
pernova remnant that is an essential element of all the models. It is however 
unclear whether the supernova cavity is required to reproduce the observed 
morphology. Another co mplication is th e non-relativistic wind of the su- 
percritical disc of SS433 ( JFabrikal . 12004 ). Its power is about an order of 
magnitude less than that of jets, but still sufficient to provide the energy of 
an average supernova explosion of ~ 10 51 erg in ~ 10 5 yr period. Though it is 
natural to expect the compact object in SS433 to be formed in a supernova 
explosion, the significance of its contribution to the power of the nebula is 
uncertain. 

Here we make some estimates of the properties that a nebula produced by 
an SS433-like object should have when evolving in a const ant- density envi- 
ronment. We consider two main energy sources: mildly relativistic jets with 
the kinetic luminosity of Lj = 2 x 10 39 erg s _1 and a photoionizing anisotropic 
EUV/X-ray radiation source with a polar-angle-dependent SED calculated 
according to lAbolmasov et al.l ( 120091 ) . There are also energy sources like the 
supernova explosion that was likely to precede the for mation of the c ompact 
accretor and a slow wind similar to that of SS433 ( jFabrikal |2004j). Both 
processes are expected to be much less energetic than the jet activity and 
photoionizing radiation of the object. 

In the next section we consider the basic properties of the HH-region 
produced by the photoionizing radiation of the central source. In section [3] 
we present several quasi-2D photoionization models. In section H] we consider 
the evolution of a jet-blown (beambag) nebula. We compare the results with 



the existing properties of ULXNe and related objects (including W50 and 
the putative relic ULX bubble in IC10) and discuss our results in section |3J 

2. Stromgren Zone 

Supercritical accretors are expected to be anisotropic sources due to ge- 
ometrical and relativistic beaming effects. Both X-ray radiation and some 
part of the EUV photons (that are most relevant for nebular physics) escape 
only in the certain range of polar angles forming a complex non-spherical 
HH-region. We introduce here a generalized version of the Stromgren radius: 

/ \ 1/3 

R s =(J- 5 ) ~ 260 n - 2/3 ^ 3 pc (1) 

Here no i s hydrogen density, a is h ydrog en recombination coefficient (see 
for example lOsterbrock and Ferlandl (120061 )). s = dS/dVt is the number of 
hydrogen-ionizing quanta emitted in a unit solid angle. S50 is s in 10 50 s _1 
units. For an isotropic source, R s is simply the radius of its Stromgren sphere. 
A strongly anisotropic source produces a bilobial or elongated nebula like the 
one shown in figure [D 

A quasi-2D sectoral approximation is more than sufficient for order-of- 
magnitude estimates. Deviations from this approximation are mainly due to 
re-emission of absorbed ionizing quanta and are therefore ~ {a a — 01 b)/® a ~ 
0.1, where a a m are recombination coefficien ts for the A and B cases, cor- 



respondingly (lOsterbrock and Ferlandl . 120061 ) . Deviations from the sectoral 



approximation are expected to be of the order 10%. As long as these effects 
are not very important one may consider the nebula extended by Rs{9) in 
any direction characterised by polar angle 9. In the next section, we calculate 
the structure of the nebula in more detail but still in a sectoral approxima- 
tion: at any polar angle 9, the radial structure is calculated as for a spherical 
nebula using our supercritical accretion disc model spectrum for appropriate 
inclination angle i = 9. 

A Stromgren z one has a characteristic expa nsion time equal to the re- 



combination time (JDopita and Sutherland! . 120031 ): 



t s = l/an e ~ 1.2 x itfTl^n^yi (2) 

-1/0 

The expansion proceeds roughly according to the law R(t) oc (1 — exp {—t/t s }) . 
Internal pressure leads to further expansion of the Hll-region but the expan- 




bow shocks 



Figure 1: Principal scheme of a hybrid nebula produced by a supercritical accretor. Pro- 
truding parts of the nebula are produced by the anisotropic radiation field. 



sion velocities are close to the thermal and are therefore overridden by the 
expanding jet-blown bubbles. We discuss this in more detail in section [H 



We use our model of an outflow-dominated supercritical disc (lAbolmasov et al. 



20091 ) to calculate the inclination-dependent spectral energy distribution. It 
is dominated by X-rays and extreme ultraviolet (EUV) radiation inside the 
funnel (for polar angles less than the funnel half-opening angle, Of ~ 20°) 
and by UV and optical at larger angles. Extended low-inclination parts of 
the nebula receive a broad-band ionizing spectrum rich in soft X-rays and in 
EUV, similar to that of active galactic nuclei (AGN). 

X-ray ionized nebulae (XINe) differ from ordinary nebulae powered by 
EUV radiation in several ways. One is that XINe do not have any sharp 
outer boundary but the boundary is determined by either absorption of X- 
rays or the equilibrium temperature. The outer radius of a XIN may be 
estimated as the radius of penetration of standard or soft X-rays bearing 
most of the energy. 

Rx ^ — ^y ( 3 ) 

where cr(E) is the absorption cross-section at energy E ~ T ~ 1 keV below 
which mos t of th e quanta are emitted. Using the simple approximation of 



Hayakawal (119731 ). one arrives to the following estimate: 



Rx ^ l.GTf^kpc (4) 

T\ here is the characteristic temperature of the energy source in keV. The 
electron temp erature fades away in XI Ne in a power-law fashion (roughly as 
T oc r~ 1//2 , see iRappaport et al.l (1994J)) therefore the definition of the outer 



radius of a nebula powered by X-rays is largely uncertain. 

3. Quasi-2D Modelling 

In order to understand the basic properties of the photoionized nebu- 
lae we suggested that the Stromgren zone is established in a quasi-2D way: 
the structure of the Hll-region in every given direction is determined by 
the generalized quanta production rate only (that is a function of the po- 
lar angle). For every polar angle value, a spherical model was calculated. 
A set of 50 spherical models for polar angles between and n/2 was cal- 
culated. Three series of spherical Cloudy (version 08.00) models were com- 
puted: one for a low ISM density (1cm -3 ) and solar metallicity, one for an 
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Table 1: Integral luminosities (10 37 ergs l units) of some selected lines for the three 
CL UD Y models described in section [3] 



line ID 


1cm 3 , solar metal- 


100cm 3 , solar 


100cm 3 , 0.1 solar 




licity 


metallicity 


metallicity 


H/3 


100 


106 


110 


[0 III]A5007 


990 


1260 


320 


[0 II]A3727 


270 


154 


60 


[0 I]A6300 


44 


39 


13 


He IIA4686 


38 


21 


20 


[S II]A6717 


28 


30 


13 


Fe VIIA6087 


5.1 xlO -3 


8.3xl0 -3 


1.8xl0 -3 


He IA4471 


8.7xl0 -3 


7.1 xlO -3 


7.4xl0 -3 


[0 III]A4363 


0.044 


0.066 


0.037 


[Fe III]A4659 


4.3 xlO -3 


2.2xl0 -3 


4.2xl0 -4 



HH-region density of 100cm -3 and solar metallicity, and one for the high 
density of 100cm -3 and metallicity ten times less then solar. The inner ra- 
dius is set to lOpc in order to reproduce the compact inner cavern created 
by shocks. T he input spectrum was calculated using the supercritical fun- 
nel model bv lAbolmasov et all fJ2009h with m = 100, 9 f = 0.4rad ~ 23°, 
r in = l/6m = 1 /600 and self- irrad i ation accounted for by two iterations (see 
original work by lAbolmasov et al.l ( 120091 ) for details). 

In figures [2rr|U we show the two-dimensional emissivity maps (cross- 
sections) of the model nebulae using combined sets of Cloudy simulations. 
The scale everywhere is linear. As may be seen in the figures, a nebula usu- 
ally has two characteristic radii, one corresponding to the Stromgren zone 
of the disc wind photosphere and the other to the nebula produced by the 
funnel interior. This structure is well established in recombination lines (see 
figure |2]) . The emissivity of a collisionally excited line is high in a narrow 
layer where the relevant ionic species is abundant and fades smoothly away 
outside the Stromgren radius where ionization by X-rays is important (see 
figure E]). Other lines such as He IA4471 and coronal emissions are sensitive 
to additional heating (see figured]). 

The luminosities of selected emission lines calculated with the three mod- 
els are given in table [TJ Note the HeIIA4686 to H/3 flux ratios about 0.2^-0.3 
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Figure 2: Model emissivities of recombination lines: H/3 (upper row) and He IIA4686 
(lower row). From left to right the pictures correspond to the low-density, high-density 
and low-metallicity models. 
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Figure 3: Model cmissivities of oxygen lines excited by collisions (from top to bottom): 
[O IJA6300, [O II]A3727 and [O III]A5007. The models are in the same order. 
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Figure 4: Model emissivities of emission lines sensitive to heating: He IA4471 (upper 
row), [O IIIJA4363 (second row) and "coronal" Fe VIIA6078 (bottom). The models are in 
the same order. 
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and very high [OIII]A5007 to H/3 flux ratios around 10 (for higher metallic- 
ity). These values together with rather bright collisionally-excited lines of 
low ionization potentials like [O I]A6300 and [S II]A6717 are consistent with 
the spectra of ULX nebulae proposed as HH-regions ex cited by the hard 



centra l radiation so urce such as the nebula of M1 01X98 (lAbolmasov et al 



2007bl ) and MH11 (lAbolmasov and Moiseevl . 120081 ). The properties of some 



individual ULX nebulae will be discussed in section EJ It should be also 
noted that the real nebular spectra are usually obtained with either a slit or 
a finite-size aperture that distorts the observed spectrum. Luminosities of 
low-ionization emissions may be underestimated by observational data. 

4. Bubble Nebulae 

4-1. Jet Propagation Law 

The expected evolution of t he nebula produced by relat ivistic jets is sim- 



ilar to the scenario proposed by lBegelman and Cioffil (Il989f ) for AGN. Unlike 



the intergalactic and galactic media in the case of radiogalaxies, the inter- 
stellar medium density is more or less constant that prevents ULXs from 
forming extended shock-powered lobes. 

The jet propagation is governed by the momentum balance at the jet 
head. In the reference system of the bow shock, the ram pressure of the 
unshocked material plus thermal pressure of the medium should be balanced 
by the momentum flux in the jet equal to ^(7 — l)Lj/2cflj where flj is the 
solid angle of a single jet. If jets are precessing or jittering, Qj effectively will 
be larger. Lj is the luminosity of the jet pair. 

We neglect the pressure of the unshocked jet material. In the case of 
SS433, the gas pressure is about eight orders of magnitude less than the 
bulk-motion momentum flux in the region where moving emission lines are 
formed. If this ratio holds, contribution of the jet gas pressure is always 
negligible. This statement is even stronger if the gas is a subject to radiative 
or adiabatic cooling. 

d_A 2 1 /3(7-l)£j ~ 

dtj ~~ 2cQjZ 2 y } 

The solution may be found analytically. There is a characteristic time tj of 
jet deceleration, and a maximal length of the jet Rj that can not be exceeded. 
Until Rj is reached, jet head moves according to the law: 
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C = vM2-r), (6) 

where ( and r are correspondingly the dimensionless jet length and age, 
( = z/Rj and r = t/tj. Final position of the jet head in this approximation 
equals: 




Here, L39, uq and T4 are the jet kinetic luminosity in 10 39 ergs _1 , inter- 
stellar medium hydrogen density in cm -3 and its temperature in 10 4 K (note 
that the gas is ionized and heated by the radiation of the central source). 
The jet head position is stationary for t > tj because the interstellar pressure 
P becomes significant. 

Characteristic time scale tj of jet deceleration equals: 

LH 2 rio 1/2 T 4 Myr (8) 



The jet bubble size is the size of a pressure-driven bubble fille d with hot 




gas an d expandi ng according t o the well- known law establi shed by lAvedisova 



(119721 ) (see also ICastor et al.l (119751 ) and lLozinskaial (119861 )): 



R JB ~ 100L^ 5 72 1/5 te /5 pc, (9) 

where t 6 is time in million years (Myr). Radius of a wind-blown bubble 
may be estimated in the same way. In the case of SS433, the power of the 
isotropic wind is about one order of magnitude lower. 

R WB ~ 70Li /5 38 n - 1/5 ^ /5 pc (10) 

The characteristic radii as functions of time for fixed ISM density and 
other parameters are shown in figure El 

4-2. Line Luminosities 

Below we use the luminosities of Balmer lines as characteristic. Two en- 
ergy sources are expected to produce strong nebular lines: photoionization by 
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the central X-ray/EUV source and radiative shock waves. One may estimate 
the contribution of the shock-powered part of the nebula as: 

Ls(HP)=(3.2QxlO-*V 2 -°- 59 + 

+4.32 x 10- 3 V 2 -°- 72 Q(V 2 - 1.5)) Lj, [ii) 

where Lj is the mechanical luminosity of the wind or jets, V 2 is the velocity 
(in 100 kms -1 units) of the shock actually producing the optical emission. 
9 here is the Heaviside function, equal to if its argument is less than and 
1 otherwise. It is used to swit ch on the shock wave precurs or for velocities 



higher that about 150 kms 1 ( iDopita and Sutherland! . Il996l ). 



The H/3 luminosity of the Hll-region may be estimated as: 

Mtf/^SxHFQsoergs- 1 (12) 

Q50 here is the total production rate of hydrogen-ionizing quanta in 10 50 s _1 
units. An EUV source with a luminosity Leuv ~ lO^ergs" 1 provides Q ~ 
4 x 10 50 s _1 . 

Several important details should be mentioned concerning emission line 
luminosities: (i) shock waves should become radiative in order to produce 
optical nebular emission; (ii) shock emission disappears if the velocity is very 
low, Vs < 20km s" 1 ; (Hi) a shock wave leaves a cavern of rarefied gas practi- 
cally transparent to ionizing radiation. Generally, photoionization produces 
brighter lines but the shape of the nebula is determined by shock waves. 

We assume that the primary source of mechanical energy is a pair of jets 
with a total power Lj = 2 x 10 39 ergs _1 in order to simulate an SS433-like 
system. Below we will trace the evolution of such a system for several million 
years. All the jet power is supposed to be distributed over a given solid angle 
flj due to precession or jitter; we assume flj = 0.038 sr. If one considers 
an SS433-like object with a stationary jet (Qj ~ 10 -3 -f- 10~ 4 ), Rj becomes 
much higher, about l-i-3 kpc. As far as I am concerned, structures of this 
kind (kiloparsec-scale double-lobe nebulae) have not ever been observed. ISM 
density is likely to vary significantly on these spatial scales. In a spiral galaxy, 
the most likely fate of a kiloparsec-scale jet is leaving the galactic disc and 
producing a fountain. 

Characteristic time scales are: the Stromgren zone establishment time, 
ts, the jet propagation time, tj, the shock cooling (to temperatures wh en the 



shocked gas becomes radiative; see for example ICastor et al.l ( 119751 )) time 

^cool- 
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t s ^ lQ^ 2 n- 1/2 T 4 - 3/4 Myr (13) 

% r l/2-l/2 




L^< 1/Z T 4 Myr (14) 



t cooi ~ 3 x 10 4 L^ 2 n 1/2 yr (15) 

In figures [MS we illustrate the dependence of various parameters of the 
nebula on time and ISM density. Generally, the nebula around a supercritical 
accretor is powered by photoionization until shock waves become radiative. 
Its properties are expected to be similar to the model nebulae described in 
the previous section, but the hot cavern in the center expands with time. 

A double bow shock nebula may exist even at longer times but usually jet 
bubble expansion velocities are higher and the nebula appears to be either a 
bubble or a multi-bubble shock-powered shell. Bubble coalescence is a poorly 
constrained process. Under almost any possible conditions, expansion of the 
two jet-powered bubbles is fast enough to produce a single cocoon at the very 
moment of jet launching. 

4-3. Double Bow Shock and Bubble Stages 

The Hll-region is torn from inside by the expanding jets and jet bubbles. 
I consider the nebula in the double bow shock dominated stage if the jet head 
propagation velocity is higher than the jet bubble expansion velocity. The 
bubble expansion is expected to be highly anisotropic in the "bow shock" 
stage but the details are uncertain and depend on the gas flows inside the 
bubble. I would also expect the morphology of the shell to depend on the 
ratio of the two velocities. I assume the nebula is in the "double bow shock" 
stage if the jet propagation velocity is higher than bubble expansion velocity, 
and in the "single bubble" stage otherwise. 

Generally, "double bow shock" nebulae exist in rarefied medium indicat- 
ing that there should be a correlation between the elongated or bilobial form 
of the nebula and the low density of the ambient medium. ULXs located in- 
side star-forming regions or molecular clouds are expected to evolve rapidly 
toward the quasi-isotropic single bubble stage. 
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Figure 5: Emission regimes of a nebula powered by a supercritical accretor. Shock waves 
are non-radiative in shaded regions (in the darkest area, both bow shocks and bubbles are 
non-radiative). The straight line with ticks demarcates the regions where photoionization 
and shock waves are expected to dominate the Balmer line emission, respectively. At the 
dotted line, the sizes of the jet-blown bubble and the Stromgren zone are approximately 
equal. 



15 



100.0 r 



l-O 

I 



U 



10.0 r 



1.0 r 



0.1 




0.01 



0.10 
T, Myr 



1.00 



Figure 6: Various characteristic radii as functions of system age and ISM density. Rs, 
Rjb and z are correspondingly the Stromgren radius, jet-blown bubble radius and jet 
length. See text for details. 
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Figure 7: Jet head and jet bubble expansion velocities as functions of the age and ISM 
density. The dotted line shows the boundary between the double bow shock case (jet prop- 
agation velocity is higher than the bubble expansion velocity) and the bubble/multibubblc 
stage. The dashed line demarcates the regions of fast (with precursors) and slow shock 
waves. In the shaded region, shocks are non-radiative (cooling time is comparable to or 
higher than the age of the system) . 
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5. Discussion 



5.1. Comparison with Observations 

The large diversity of the observational properties of ULX nebulae may be 
explained (at least qualitatively) in the framework of a hybrid model taking 
into account both EUV radiation and shock waves. Let us see how individual 
objects fit into this unified scheme. 

MF16 is a well-known but outstanding example. It is the most com- 
pact of all the known ULXNe and has the highest expansion rate measured 
(V ~ 200km s" 1 ). It was shown to have an elongated bilobial morphology 
with a fainter extended halo su rrounding a prominent shell with a radius 



of 10 -r 20pc (jBlair et al.l . 120011 ). All these parameters are consistent with 
the position around Uq ~ 5 cm -3 , t ~ 30 000yr in our diagrams. The ex- 
tended halo probably corresponds to the outer part of the Hll-region that 
is expected to be larger than the jet bubb l e rad ius at this age. The ISM 
density was estimated in lAbolmasov et al.l (120081 ) and equals no — 6 cm -3 . 
The uniqueness of MF16 is probably the result of its evolutionary stage that 
lasts less then one tenth of th e expected lifetime of the source. 

Several objects selected by lPakull and Mirionil (j2003h are large-scale bub- 
bles with dynamical ages ~ lMyr and quasi-spherical morphology probably 
affected by underlying density gradients. It is likely that their expansion 
started in a rather dense environment of the host star-forming region that 
implies an effective density value of no ~ 10 -j- 100 cm -3 . 

Predicted lifetimes for nu clear-timescale ac c retion in massive X-ray bi- 



naries are of the order lMyr (IRappaport et al.l . 120051 ) . which is comparable 



to the dynamical ages of ULX bubbles. It is possible that the bubble sur- 
vives the further evolution of the object that is expected to evolve towar d 
a binary system similar to Cyg X-3 or IC10 X-l (JBauer and Brandtl . 12004 ). 
consisting of a black hole and a Wolf-Rayet donor. In this scope, th e giant 
synchrotron superbubble around IC10 X-l (JLozinskaya and Moiseevl . 120071 ) 
is the best candidate for a relic ULX bubble. Its ki nematical p r opert ies are 
close to those of the largest ULX bubble nebulae (lAbolmasovl . 120081 ) . The 
concurrent explanation of the bubble in IC10 as a hypernova remnant can 
not account for this similarity. 

HoII X-l does not show a ny signatures of supers onic motions or contam- 



ination with heavy elements (JLehmann et all 120051 ). The nebula appears to 



be large r then the Stromgren radius inferred from recombinati on line lumi- 
nosities (lAbolmasov and Moiseevl . 120081 ; lAbolmasov et al.l . l2008l ) that may be 



explained by a non-radiative central cavern created by the preceding super- 
nova explosion and/or wind/jet activity. This is consistent with the observed 
non-thermal radiospectrum of the nebula. 

W50 is also quite easy to p lace on the diagram. The undisturbed ISM 
density ~ 2 cm -3 is known from lLockman et al.l (J20071 ) so the nebula may be 
considered a prototype of an evolved ULXN in a rarefied environment. Its 
properties are consistent with an age about several xl0 5 yr and ISM density 
~ 2 cm -3 . These estimates imply that the bubble expansion velocity for 
W50 should be Vjp ~ 50 ± 100km s" 1 consistent with veloc i ty est imates by 



Zealev et al.l (119801 ): iBoumis et al.l (120071 1: lAbolmasov et al.l (12010h . 



5.2. Jet-Powered Nebulae 

The ambiguity of the energy sources of ULXNe bears some similarity 
with the question about the power sources of Sey fert Narrow-Line Regions 
(NLRs). This similarity is most striking for MF16 (lAbolmasov et al.l . l2007af ) 
where the integral spectrum itself is similar to that of a Seyfert NLR. 

Another parallel with AGN is in the pair of relativistic jets that in the case 
of SS433-like object s are likely to produce som ething similar to the cocoons 
that were predicted (jBegelman and Cioffil . |l989j) but are not observed around 
AGN mainly due to strong gradients in ISM and IGM density. 

Little can be said for certain about the structure of a nebula produced 
by relativistic jets. However, W50, ULXNe and similar objects are prospec- 
tive targ ets for sophisticated hydrod yna mical modeling . For example, in the 
works by lVelazquez and Ragal (120001 ) and lZavala et al.l (120081 ). hydrodynam- 
ical simulations allowed to qualitatively explain the basic properties of W50. 
I suggest that running simulations at slightly different parameters would al- 
low to better understand the diversity of ULXNe that sometimes evolve into 
quasi-spherical bubbles, but under different conditions may retain a resem- 
blance to W50 instead (as in the case of MF16). 



6. Conclusions 

In this article, I have shown that most of the observational properties of 
ULX nebulae may be explained by a single unified scheme. Only two free 
parameters were varied (ISM density and the age of the source), though some 
divergence in the parameters of the binary system itself is expected. 

One of the main results is the bipolar morphology in some "hot" lines 
such as He IIA4686 and [O III] emissions. The "ionization cone" structures 
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detected by iRoberts et al.l (120031 ) and the elongated high-excitation nebulae 
around M101P098 and HoIX X-l are naturally explained in the framework 
of our model. Bilobial and elongated shock-powered nebulae are expected 
to appear in a low-density environment due to the different dependences of 
the jet propagation and bubble expansion velocities on the ambient density. 
The morphology of ULX shells is predicted to evolve through a "double bow 
shock" bipolar structure toward a quasi-isotropic single bubble stage. 

Supercritical accretors are likely to exist for about lMyr(of the order of 
the nuclear timescale of the donor star). Then the source is expected to 
evolve into a WR+BH binary, possibly surrounded by a relic ULX bubble. 
We propose that the synchrotron bubble in IC 10 is a relic ULX bubble rather 
than a hypernova remnant. 
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